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1. Introduction 

In the last years, new measurements of heavy-quarkonium decay observables, mainly com- 
ing from BABAR, BELLE, BES, CLEO and the Fermilab experiments have improved our 
knowledge of inclusive, electromagnetic and several exclusive decay channels as well as of 
several electromagnetic and hadronic transition amplitudes. The new data and improved 
error analyses of the several correlated measurements have not only led to a sizeable reduc- 
tion of the uncertainties but also, in some cases, to significant shifts in the central values 
||l|. Such data call for comparable accuracies in the theoretical determinations. 

The main mechanism of quarkonium decay into light particles is quark-antiquark an- 
nihilation. It happens at a scale which is twice the heavy-quark mass M. Since this scale 
is perturbative, quark-antiquark annihilation may be described by an expansion in the 



strong-coupling constant Og- Experimentally this fact is reflected into the narrow widths 
of the quarkonia below the open flavor threshold. The bound-state dynamics, instead, is 
characterized by scales that are smaller than M. While we may not necessarily rely on an 
expansion in Os to describe it, we may take advantage of the non-relativistic nature of the 
bound state and expand in the relative heavy-quark velocity v. Decay-width formulas may 
be organized in a double expansion in Og, calculated at a large scale of the order of M, and 
in V. In the bottomonium system, typical reference values are as{Mh) ~ 0.2, v'^ ~ 0.1 and 
in the charmonium one, as{Mc) ~ 0.35, f^ ~ 0.3. We will call relativistic all corrections of 
order v or smaller with respect to the leading decay-width expression. 

S-wave quarkonium decays into lepton pairs are known within a few percent uncer- 
tainty. Theoretical accuracies of about 5% in the charmonium case and of about 1% in 
the bottomonium one require the calculation of 0(f^, Og v^jOg) corrections. Two photon 
decays are not so well known. The rjc -^ 77 width is known with an accuracy of 35%, the 
XcO -^ 77 width with an accuracy of 20% and the Xc2 — > 77 width with an accuracy of 10% 
Q. Nevertheless, in the P-wave case, the improvement has been dramatic over the last few 
years: the errors quoted in the 2000 edition of the Review of Particle Physics were 70% for 
the XcO -^ 77 case and 30% for the Xc2 -^ 77 case [Q . Theoretical accuracies matching the 
experimental ones require the calculation of 0{v'^,as) corrections. 

In this work, we consider relativistic corrections of order f ^ and order v'^ to electro- 
magnetic decays of S and P wave respectively. Since the leading-order S-wave decay width 
is proportional to the square of the wave function in the origin, it is of order v^. Since 
the leading order P-wave decay width is proportional to the square of the wave-function 
derivative in the origin, it is of order v^. Therefore, both corrections of order v^ to S-wave 
decays and of order v'^ to P-wave decays provide decay widths at order v'^ . 

In the S-wave case, corrections of order v'^ and f^ were first considered in [^ and 
in respectively. We agree with their results if we use their power counting, but we find 
three new contributions using our different counting. Moreover, we resolve an ambiguity 
in the matching coefficients of the operators contributing at order u^. In the P-wave case, 
corrections of order v'^ were first calculated in P]. Our results disagree with those in H. 
We find that we need two operators less to describe the decay widths at that order. 

This work is partially based on [0]. We refer to it for some detailed derivations. 

The paper is organized in the following way. In Sec. |2|, we set up the formalism, 
discuss the power counting and introduce our basis of operators. In Sec. |^, we perform 
the matching. We also discuss the cancellation of the infrared singularities in the matching 
of the octet operators. In Sec. Q, we conclude by discussing applications and further 
developments of this work. In Appendix ^ and ^ we list all the operators and the matching 
coefficients that have been employed throughout the paper. 

2. Decay widths in NRQCD 

2.1 NRQCD 

In the effective field theory (EFT) language of non-relativistic QCD (NRQCD), the annihi- 
lation of quarkonium is described by four-fermion operators Q. The decay width factorizes 



in a high-energy contribution encoded in the imaginary parts of the four-fermion matching 
coefficients and a low-energy contribution encoded in the matrix elements of the four- 
fermion operators evaluated on the heavy-quarkonium state. The NRQCD factorization 
formula for the electromagnetic decay width of a quarkonium H is 

TiH -. em) = 2^ ^(^l^^&Ll^)- (2-1) 

n 

In this work, H -^ em stands either for H ^ jj or for H — > e^e^. \H) is a dimension 
—3/2 normalized eigenstate of the NRQCD Hamiltonian with the quantum numbers of 

(n) 

the quarkonium H. O^^^^^ is a generic four-fermion operator of dimension d„. The "em" 
subscript means that we have singled out the electromagnetic component by projecting onto 
the QCD vacuum state |0). The general form of 0^_^^^^ is '0^(- • • )x |0) (0| xH' ' ' )V') were ip 
is the Pauli spinor that annihilates a quark and x is the one that creates an antiquark. The 
operators (• • • ) may transform as singlets or octets under SU(3) gauge transformations. A 
list of four-fermion operators is provided in Appendix ^. The coefficients Cem encode the 
high-energy contributions to the electromagnetic annihilation processes, which have been 
integrated out from the EFT. At variance with Q, we will not label differently coefficients 
that stem from decays into 77 or e~^e~. The coefficients are calculated in perturbation 
theory by matching Green functions or physical amplitudes in QCD and NRQCD. It is the 

(n) 

purpose of this work to calculate Imcem for operators up to dimension 10. 

2.2 Power Counting 

In the factorization formula (|2.l|), the matching coefficients Cem are series in a^ while the 



matrix elements {H\0^^^^\H) are series in v and are, in general, non-perturbative objects. 
It is not possible to attribute a definite power counting to the matrix elements because 
of the several contributing energy scales, as it is typical in a non-relativistic system: Mv, 
Mv'^, ... . Whatever power counting one assumes, as long as v <C 1, operators of higher 
dimensionality are suppressed by powers of v. 

In the following, we will assume Mv to be of the same order as the typical hadronic 
scale Aqcd and adopt the following counting. Matrix elements of the type {H'\0\H), where 
0\H) and \H') have the same quantum numbers and color transformation properties in the 
dominant Fock state, scale (at leading order) like {Mv)'^~^, d being the dimension of the 
operator O. liO\H) and \H') do not have the dominant Fock state with the same quantum 
numbers, then the matrix element singles out a component of the quarkonium Fock state 
that is suppressed. The amount of suppression depends on the power counting and on the 
quantum numbers. Let us consider a component in \H') that shows up as a first-order 
correction to a singlet quarkonium state with orbital angular momentum L, spin S and 
total angular momentum J. The correction is induced by an operator T of dimension dj^: 



S',L',J' 



l2S'+lTl I ^ T7|25+lr \ 

E-E' 



E 1^^'"^^'.')^ — T:t;' "^^ (2.2) 



This gives a nonvanishing contribution to the matrix element {H'\0\H) if for some S", L' 
and J' the matrix element in ( p.2| ) does not vanish and if there is a Fock state of 0\H) with 
the same quantum numbers as p ^^Lj,). In case this is the dominant Fock state of 0\H), 
the suppression factor with respect to the naive scaling, {MvY~^ , is Mv'^^~'^ /{E — E'). 

As an example, we consider the dimension 5 operator T = tpH t"" gA"" ■ dip + c.c, where 
t°' are the color matrices in the fundamental representation and c.c. stands for "charge 
conjugated" {ip -^ ic^X*-, *° ^)! ^ ~(*°^^)"^)- Inserted in (|2l^ , it projects onto an 
octet state with quantum numbers L' = L it 1 and S' = S. Color octet quark-antiquark 
states may appear in the quarkonium spectrum combined with gluons to form hybrids. 
Gluonic excitations of this kind are expected to develop a mass gap of order Aqcd, hence 
E — E' ^ Aqcd ~ Mv. As a consequence, in a quarkonium state the octet component with 
quantum numbers S and L it 1 is suppressed by v with respect to the singlet component 
with quantum numbers S and L. Note that to generate a D-wave component {L = 2) 
from a S-wave state we need two operator insertions and second-order perturbation theory. 
Hence, this is suppressed by v"^ with respect to the S-wave component. 

A similar counting holds if we consider the dimension 5 operator J^ = ip^ t"" gB"^ ■ aip + 
c.c. Inserted in ( |2.2| ), it projects onto an octet state with quantum numbers L' = L 
and S' = S lb 1. Thus, also the octet component with quantum numbers S" it 1 and L is 
suppressed by v with respect to the singlet component with quantum numbers S and L. 

The power counting adopted here is the most conservative one in the framework of 
NRQCD. This is the reason for our choice. The price we pay is that some observables will 
depend on more matrix elements than they would in a different counting. An example of 
alternative power counting, which seems better suited for the situation Mv^ ~ Aqcd, is 
provided by Ref . [Q] . For a critical review and a discussion on the different power countings 
we refer to |^] and references therein. 

2.3 Four-fermion operators 

The four-fermion sector of the NRQCD Lagrangian contains all four-fermion operators 
invariant under gauge transformations, rotations, translations, charge conjugation, parity 
and time inversion. They may be classified according to their dimensionality and color 
content. Moreover, it is useful to decompose them in terms of irreducible spherical tensors. 
Some of the operators are redundant because they may be expressed in terms of others by 
means of field redefinitions. Finally, the power counting introduced in the previous section 
sets the relevance of the different operators in the decay width formulas. We will address 
all these issues in the following. 

2.3.1 Operators from dimension 6 to 10 

The NRQCD operators may be organized according to their dimensionality. Four-fermion 
operators of dimension 6 were considered in Q. For dimensional reasons, they cannot 
depend on the gluon fields. The only allowed color structures are \c ® Ic and t"" ® t"". 
However, since we are considering electromagnetic decays whose final state is the QCD 
vacuum, the t" ® f^ structure is forbidden. 



Parity conservation forbids four-fermion operators of dimension 7. Four-fermion op- 
erators of dimension 8 may be built with two covariant derivatives [Q] , or with a chromo- 
magnetic field, like, for instance, 

V^t^fi-axlOXOlxV + H.c, (2.3) 

where H.c. stands for Hermitian conjugated. The operator describes a (spin-flipping) octet 
to singlet QQ pair transitions. The same operator with a chromoelectric field instead of 
the chromomagnetic one is forbidden by parity conservation. 

Operators of dimension 9 may involve a covariant derivative and a chromoelectric field, 
like, for instance, 

V^t^D •ax|0)(0|x^5^-<7V^ + H.c.. (2.4) 

The same operator with a chromomagnetic field instead of the chromoelectric one is for- 
bidden by parity conservation. 

Dimension 10 operators may involve either four covariant derivatives, or two covariant 
derivatives and a chromomagnetic field, like, for instance, 

V^t5^-'Dx|0)(0|x^'D -aV^ + H.c., (2.5) 

or two gluon fields. 

2.3.2 Irreducible spherical tensors 

All the tensorial structures that are consistent with the discrete symmetries and rotational 
invariance are allowed. So, for instance, besides Eq. (|2.3| ) also the following operator is 
possible 

V'^ff^x|0)(0|x^a^ + H.c.. (2.6) 

It is useful to decompose Cartesian tensors into irreducible spherical tensors. This 
allows a classification of the operators in terms of angular-momentum quantum numbers: 
~^^Lj, where S is the spin, L the orbital and J the total angular momentum. 

According to3(8)3 = 503©1, a tensor made of two vectors A and B may be 
decomposed as 

A-BJ = A('B.) + i;^i:ii^ + ^i.B, (2.7) 

2 3 ^ ^ 



where 



A^'B^^ = 1 ^A . i?, (2.8) 



is a symmetric and traceless tensor belonging to the representation L = 2 of the rotational 
group, and the last two terms of Eq. ( |2.7] ) are tensors transforming respectively as L = 1 
and L = spherical harmonics. So, for instance, besides Eq. ( |2.4D also the following 
operator has to be considered 

^t'D (VJ')x|0)(0|x^5-E(VJ'V + H.c. (2.9) 



Note that a symmetric tensor S^^ may be decomposed, according to 6 = 5 1, into 

S'i = S^'i^ + -^5^^ (2.10) 

where S^*-^' is a traceless symmetric tensor. 

Therefore, from (5©1)(8'3 = 705©3©3, a tensor made of a symmetric tensor S"^^ 
and a vector A may be decomposed into 

Jklxjm <_ ^jklxi'm I JmpQ(np)An\ 

were 5^^*-^^^ stands for a tensor symmetric with respect to the indices i, j and k, and 
such that all partial traces (obtained contracting two of the three indices) are equal to zero: 

g{(ij)j^k) _ £ ( s^'^^^A'' + 5'(*'^)yl-' + S^^'^^A^^ — (5^^^'^ + S^'^6^^ + 5^''6^^') S^'^^'^A"^. 

(2.12) 
The second term in Eq. ( p.ll| ) transforms like a L = 2 spherical harmonics. Note that 
^inpg{mn)j^p _j_ ^imp g{np) jiji j 2 jg symmetric in the indices I and m and traceless. The third 
and fourth terms in Eq. ( 2.1l| ) transform like vectors. 



As an application, we may use Eq. ( 2.11| ) to decompose the tensor 



4°') [-in v^") <'■"' 

into irreducible spherical tensors. If we treat the covariant derivatives as ordinary ones, 
the tensor, being completely symmetric, has 10 independent components. According to 
10 = 7 © 3, it may be decomposed into a L = 3 tensor 

Tgi^ (4)^»-«" (41?') - I (4) «"^'" (-5«)^ (2.14) 

where A^^J^^ = A^5^^/2 + A^^'' /2 - A^5'^ /^, and a tensor that transforms like a L = 1 
representation of the rotational group: 

From these we may construct the four-fermion operators listed in Eq. ( |A.9D : 

VU'Po) = ^V^T;jyx|0)(0|xt^ (-^'^•<?) V' + H.c, (2.16) 

VemeP2) = ^^t M'^J^' -^-^T\ll\ a\|0)(0|xt(-^^(V))^ + H.c.,(2.17) 

VemeP2,'F2) = ^V^xg^ x|0) (Olx^ (^_1^(V^) j ^ + H.C.. (2.18) 



2.3.3 Field redefinitions 

Operators may be redundant, in the sense that they may be traded for others by means of 
suitable field redefinitions. Let us consider the following cases. 

(i) First, we consider the field redefinitions (a is a free parameter) 



Tp ^ ifj + 



M5 



X^X 



M5 



D] ,xlO)(0|x^ 



'-"d] ,v|0)(0IV'^ 



^ 



X 



(2.19) 



which induce the following transformations: 



V''^ iDq tp + xUDq X ^ i^^ iDo ip + x^ iDo x 






M5 



'?8cm( 5*0)) 



2M^ ^ 2M 



(2.20) 



i^^^i^-X^^X + 2^ {Q'en.eSo) - Q"em(^5o)) , (2.21) 



where the operator 78om(^5'o) has been defined in Eq. (| 
and Q"cmC'So) in Eq. 



and the operators Q'om(^'S'o) 
In Eq. ( 2.21 ) we have neglected operators proportional to 
the center of mass momentum and octet operators, whose matrix elements are subleading 
in any power counting. Clearly, for a suitable choice of the parameter a we may elimi- 
nate the operator '78em(^'S'o) from our basis of operators in exchange for a redefinition of 



the matching coefficients of Q'em(^5'o) and Q"em(^5'o) 
field redefinition. 



hL 



{'So) 



^em(^'S'o) + 2a and 



h'lr^i^So) —f /i"jn(^S'o) — 2a. Note that the sum h'^^(}So) + h'^^(}So) is invariant under the 



(ii) Second, we consider, the field redefinitions 



^r^ + w^'i'L-' 



X 



X 









^^J,x|o)(o|x^ 
'dA,V'|o)(o|v^ 



(jV 



<y^x 



(2.22) 



where 



.(0) ^SW^ 

" ijlk o ' 

-,(1) ^ijn^klr. 



" ijlk 2 ' 

,(2) _ 6'^6^'' + 6^^6'^ 6'^ 5^^ 

-ijlk ~ 9 o • 



(2.23) 
(2.24) 

(2.25) 



They induce the following transformations: 

V't iDo ^l^ + x^iDoX^ i'^ iDo ^ + x^iDox-^ ^SemC'^i), (2.26) 

+2^ (Q'em('5i) - Q"U'Si)) , (2.27) 

f^2 jj2 jj2 f)2 

+— -g X (octet operators), (2.28) 

jj2 jj2 f)2 f)2 

2M '2M^ J^ 2M 2M^ 

+2^ {Q'em^Si,^ Di) - Q",^{^Si,^ Di)) , (2.29) 

where the operators Tg^J^i^Si) have been defined in Eq. ( [A.7D and the operators Q'cm(^'S'i)> 
Q"em(^5i), Q'cmC^Si,^Di) and Q" eraC^ Si,^ Di) in Eq. (Q. InEqs. (|3^)-(|3^ we have 
neglected operators proportional to the center of mass momentum and octet operators. 
Again, for a suitable choice of the parameter a we may eliminate the operators '^semi^^^) 
from our basis of operators in exchange for a redefinition of the matching coefficients of 
Q'emC'Si), Q"em('5i), Q'em('5i,' I?i) and Q" cm^ Si,^ D^) . Note that K^{^Si) + h'U^Si) 
and h'^^{^Si,^ Di) + h"^{^Si,^ Di) are invariant under the field redefinitions. 

In 1^], it was pointed out that the operators '78em(^'S'o) and '78om(^5'i) may be eliminated 
in favor of Q'em(^'S'o) and Q"em(^'S'o), and Q'em(^5'i) and Q"em(^5'i) through the use of the 
equations of motion. Our argument, which uses field redefinitions instead of equations of 
motion, is equivalent Q. 

(in) Finally, we may wonder if there are field redefinitions similar to ( |2.19| ) and ( |2.22| ) 
that allow to eliminate from our basis also the operators T^^m^Pj) defined in Eq. 
Indeed, such redefinitions exist and are 






^r x-^,%L^\—^^']mm'<y\—^D-\x 



(2.30) 



However, in this case, they induce the following transformations 

-2]|5T(/iz5o (V'V^ (-1^^) x) |0)(0|xta' (-^^') ^' (^-^^^ 

f)2 f)2 f)2 f)2 

^^ — ^ - y^ — y ^ ^"^ — -0 - y"^ — y - 2 — aVcrni^Pj), (2.32) 

2M 2M J 2M 2M M^ ^ ' ^ ' 

where in Eq. ( 2.32| ) we have neglected operators proportional to the center of mass momen- 
tum and octet operators. We see that the field redefinitions ( p. 30 ) may be useful to trade 



the operators %,m{^Pj) for tJ/^ id^ U^a' (-0A x) |0) (0|xV' (-0n V', but, dif- 



-ijik '^i> y^y "^ y 2 -^ y ^y 1^/ \"i^ '^ V 2 

ferently from the previous cases, they do not reduce the number of operators. We wih keep 
Tgeml^-Pj) in our basis of operators. 

2.3.4 Power counting of the four-fermion operators 



From the rules of Sec. 2.2, it follows that 



{He'+'Lj)\j^ 0,me'^'Lj)\He'+'Lj)) ~ Mv'^-^ (2.33) 

where \H(^ ^^Lj)) stands for a quarkonium state whose dominant Fock-space component 
is a QQ pair with quantum numbers S, L and J, OcmC ^^Lj) is a singlet four-fermion 
operator that acts on the QQ pair like a spin S, orbital angular momentum L and total 
angular momentum J tensor and d is the dimension of the operator. 

Concerning the power counting of the octet matrix elements, first, we consider the 



matrix elements of the dimension 8 operators defined in Eq. ( A. 6 ): 



{HeSo)\^SsemeSo)\HeSo)), (2.34) 

and 

{Hes,)\^^ssemes,)\Hesi)). (2.35) 

The operator 58cm(^5'o) destroys a singlet QQ pair with quantum number ^Sq and creates 
an octet QQ pair with quantum numbers '^Si and a gluon (and viceversa), the operator 
'58em(^'S'i) destroys a singlet QQ pair with quantum number ^5i and creates an octet 
QQ pair with quantum numbers ^Sq and a gluon (and viceversa). Following the power 
counting of Sec. 2.2 the chromomagnetic field scales like (Mv)"^, moreover the octet Fock- 
space component is suppressed by v. Hence, both matrix elements scale like Mv^. 



Equations ( A. 7] ) define octet operators of dimension 9. The operator Tg^l^{^Si) de- 
stroys a singlet QQ pair with quantum number ^Si and creates an octet QQ pair with 
quantum numbers ^Pi and a gluon (and viceversa), the operators Tg em (^-Pj) destroy a 
singlet QQ pair with quantum number ^Pj and create an octet QQ pair with quantum 
numbers ^Si and a gluon (and viceversa). Following the power counting of Sec. ^^, the 
chromoelectric field scales like (Mv)'^ and the covariant derivative like Mv, moreover, since 
the octet Fock-space component is suppressed by ?;, we have 



{HeS,)\^T,ZeSl)\HeS^)) r. Mv\ (2.36) 

and 



{H{^Pj)\^%U^Pj)\H('Pj)) ^ Mv-^ . (2.37) 

Matrix elements involving dimension 10 octet operators are negligible at order v^ . An 
example is the matrix element 

(F(3po)|^ tlj^B ■ 'Dx|0)(0|x^'d • ai:\H{^Po)) ~ Mv\ (2.38) 



We have kept ambiguous the ordering of the covariant derivatives appearing in some 
of the singlet operators of dimension 10. For instance, D ^ could stand either for ( D '^)'^ or 
D^ D^ D^ OT D^D^D^D^ov combinations of them. The ambiguity can only be resolved 
by octet operators of dimension 10, which are beyond the present accuracy. 

All operators involved in the matching are listed in Appendix ^. 

2.4 Electromagnetic decay widths 

Having assumed a power counting and having chosen a basis of operators, we are in the 
position to provide explicit factorization formulas for S- and P-wave electromagnetic decay 
widths. Up to order f^, these are 

+ ^'"^tf^'\ Hesmu's,)\Hes,)) + ^-^^^^%f^{Hes,)\Q'Les,)\Hes^)) 



(2.39) 



TiHCSo) - 77) = ^^"^;;S^'\ Heso)\ou'so)\Heso)) 

+ ^I^g-/^°) (ff(is,)|p,,(is,)|ff(i5,)) + '^"^^^:^'^°\ ff(^So)|^8e^(^5o)|g(^So)) 

+ ^^^^^^%f^{HeSo)\QU'So)\HeSo)) + ^^"''^f^'\ HeSo)\QU'So)\HeSo)), 

(2.40) 

+ ^^"';f^'\ HePo)\V.^ePo)\HePo)) + ^^"'g^'^°\ ff(^^0)l?'8en.(^P0)|g(^P0)), 



(2.41) 



TiHep,) -. 77) = ^^^^^^{HeP2)nn.eP2)\HeP2)) 






+ ^-^^^^^^^^{Hep,)\vu'P2)\HeP2)) + ^''^'iTf'''\ Hep,)\Tsu'P2)\HeP2)). 



(2.42) 
In (2.39) and (2.4C|), the first matrix element scales like v^, the second one is suppressed 



by f^, the third one by v^ and the other ones by v^. In ( ^.4lD and ( 2.42 ), the first matrix 
element scales like v^ and the other ones are suppressed by v'^. 

Some comments are in order. In Ref. j^, the matrix elements of the operators 
'58em(^'S'i)) 'Pemi^Si,^ Di), Tg^^(^Si) and 58em(^'S'o) Were not included in the expressions 



10 



e^e ) and T{H(^So) — > 77) The reason is that in the power counting 



of TiHi^Si) 

adopted there, which is different from the one discussed in Sec. 2.3.4, they were considered 

to be suppressed.^ 

In Ref. H, also the matrix element 



{HeP2)\gu'P2)\HeP2)) 



of the operator 

QLeP2) 



2^ 



--) D^'D^^ 



-1^ 
2 



^]x\o){o\x^ 



--'d(V^) ) V' + H.c. 



was considered to contribute to Eq. (2.42). In the basis of operators given in (A.^ 
be rewritten as 



(2.43) 



(2.44) 



it can 



(2.45) 



The operator 'Pcm(^-P2) contributes, indeed, to the decay width of a ^P2 quarkonium and 
its contribution has already been accounted for in Eq. ( 2.42| ). The question is whether the 
operator 7^em(^-P2,^ -^2) contributes at relative order v^ to the decay width of a quarkonium 
whose dominant Fock state has quantum numbers ^^2- In Eq. ( 2.18| ), we have seen that 



the operator acts on one of the quarkonium states with a term proportional to T--la, 



.(3) 



ijk 



where T^ .^ transforms like a spherical tensor in the L = 3 representation of the rotational 

group. The Wigner-Eckart theorem guarantees that the matrix element {Q\'T\-l,a^\^ P2) 
vanishes, because we cannot generate a L = tensor by combining a L = 3 with a L = 1 
tensor. Explicitly, this is reflected by the fact that 



(0|T 



(3)^fc|3i 
ijk 



a 



d pTv 



(y) 



p ■ a 



-p'^p'^^a^- 



'n-Lnk f \\ „k 



M)P' 



0, 



(2.46) 



where hl^p (A) is a symmetric and traceless rank-2 tensor that represents the polarization 
of the ^P2 state and the trace is meant over the Pauli matrices. 



3. Matching 

In this section, we calculate the imaginary parts of the matching coefficients that appear 



in Eqs. ( |2.39| )-(2.42). The method consists in equating (matching) the imaginary parts of 
scattering amplitudes in QCD and NRQCD along the lines of Ref. Q. 

In the QCD part of the matching, the ingoing quark and the outgoing antiquark are 
represented by the Dirac spinors u{p) and v{p) respectively, whose explicit expressions are 



u{p) 




v{p) 



lEp + M 



2 En 



p-a 



(3.1) 



where Ep = yp"^ + M^, and ^ and rj are Pauli spinors. In the NRQCD part of the matching, 
the ingoing quark and the outgoing antiquark are represented by the Pauli spinors ^ and 
r] respectively. 

^We thank G. Bodwin for communications on this point. 
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Figure 1: QCD Feynman diagrams describing the amplitude QQ -^ 77 -^ QQ at leading order. 

Since we need to match singlet and octet-singlet transition operators, we consider both 
the scattering amplitudes QQ -^ QQ, with two photons or an electron loop as intermediate 
states, and QQ g -^ QQ, with two photons or an electron loop as intermediate states and 
an external gluon in the initial state. 

In the center of mass rest frame, the energy and momentum conservation imposes the 
following kinematical constraints on the scattering QQ -^ QQ, 



\P\ = \k\, 



p + p' 



0, 



k + k' = 0, 



(3.2) 



and on the scattering QQ g — *• QQ, 



Ep + Ep 



+ \q\ = 2Ek 



p+p' + q = 0, k + k' = 0, 



(3.3) 



where p, p' are the ingoing and k, k' the outgoing quark and antiquark momenta, while 
q is the momentum of the ingoing gluon, which is on mass shell. Note that, in the non- 
relativistic expansion, the gluon momentum |g| is proportional to (three-momenta)^/M. 

The matching does not rely on any specific power counting and can be performed 
order by order in 1/M |1^. We will perform the matching up to order 1/M^, which is the 
highest power in 1/M appearing in Eqs. ( 2.39| )-( p.42 ). In practice, we expand the QCD 
amplitude with respect to all external three-momenta. In the case of the QQ g -^ QQ 
scattering, the expansion in the gluon momentum may develop infrared singularities, i.e. 
terms proportional to l/\q\- These terms cancel in the matching. We will discuss the 
cancellation in Sec. |3.5| . 

Finally we note that, since the matching does not depend on the power counting 
and the scattering amplitudes do not have a definite angular momentum, the matching 
determines many more coefficients than needed in Eqs. (|2.3£ 



(2.42) 



3.1 QQ 



77 



The matching of the QQ -^ 77 amplitude is performed by equating the sum of the imagi- 
nary parts of the QCD diagrams shown in Fig. |l| (taken by cutting the photon propagators 
according to l//c^ — > —2t: i 5{k'^)9{k^)) with the sum of all the NRQCD diagrams of the 
type shown in Fig. |2[ 
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Figure 2: Generic NRQCD four-fermion Feynman diagram. The empty box stands for one of 
the four-fermion vertices induced by the operators listed in Appendix 0, Eqs. (A^)-(A^) and 



-(A.IO) 



We find that 

(3.4) 
(3.5) 
(3.6) 
(3.7) 

(3.8) 

lmh'^^{^So) + lmh'^^i^So) = ^a'^Q^n, (3.9) 

(3.10) 
(3.11) 

(3.12) 

where a is the fine structure constant and Q the charge of the quark. The four-fermion 
operators to which the matching coefficients refer are hsted in Appendix pi 

The coefficients of the operators of order 1/M^ and 1/M^ agree with those that may 



Im/e^(i5o) = 


a Q TT, 


Im5rem(^5'o) = 


4 2^4 

--a Q vr, 


Im/en.('Po) = 


3a^Q\, 


Im/en.('P2) = 


-^a'Q\ 


Im^cm(^I)2) = 


A«^gV, 


+ lmhU'So) = 


^«^gv, 


Imgcm(^^o) = 


-7a^Q\, 


Imgem(^-P2) = 


8 2/^4 

5 


a5em('P2,'i^2) = 


20 2n4 

21 ^ 



be found in Q. Here and in the following, we also refer to [11| and references therein for 
an updated list of the imaginary parts of the four-fermion matching coefficients appearing 
at order 1/M^ and 1/M^ in the NRQCD Lagrangian. Some of them are known at next- 
to-leading and next-to-next-to-leading order. The coefficient ImhemC^^) agrees with the 
result of Ref. [|l2[ . Equation (^]^) agrees with the one found in |^ . By matching the dia- 
grams of Fig. |l|we cannot resolve lm.h'^^{^So) and Im.h'^^(^So) separately. Equation ( |3.1C| ) 
agrees with the one found in Q, whereas our evaluation of the coefficients Ira gcmi^ P2) and 
lTagemCP2,^ F2) disagrees with that one in p. 

3.2 QQ -^ e+e- 

The matching of the QQ — > e^e~ amplitude is performed by equating the imaginary part 
of the QCD diagram shown in Fig. |3| (the imaginary part of the e^e~ pair contribution to 
the photon's vacuum polarization is —a {k'^g'^'^ — k'^k'^)/3) with the sum of all the NRQCD 
diagrams of the type shown in Fig. ^ 
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Figure 3: QCD Feynman diagram describing the amplitude QQ 



QQ at leading order. 



We find that 






lmg,^i^Su^D^ 



3 ' 
9 ' 



lmh'^^{^Si) + lmh"^{^Si) = — o^QV, 



Irnh'^^i'^Su'' Di) +lmh'^^{''Si;' D 



Im/lem(^I?2) =0, 
Im/lem(^I?3) =0, 

23 
" 36 



54 



12 



'/ ^3 c 3 



a^Q^vr. 



(3.13) 
(3.14) 

(3.15) 
(3.16) 

(3.17) 

(3.18) 
(3.19) 

(3.20) 



The four-fermion operators to which the matching coefficients refer are listed in Appendix 

The coefficients of the operators of order 1/M^ and 1/M* agree with those that may 
be found in [^. Equation ( p.l6| ) agrees with the one found in ||] and Eq. ( 3.17| ) with 
the D-wave decay width calculated in [12|. Equation (3.20|) is new. Equations (|3.1^) 



and ( 3.19 ) follow from angular momentum conservation. Similarly to the case discussed 
above, we note that by matching the diagram of Fig. y we can only determine the sums 
lmh'^^{^Si) + lmh'^^{^Si) and lmh'^^{^Si,^ Di) + lmh'^^{^Si,^ Di) but not the individual 
matching coefficients. 



3.3 QQg 



77 



The matching of the QQg — > 77 amplitude is performed by equating the sum of the 
imaginary parts of the QCD diagrams shown in Fig. with the sum of all the NRQCD 
diagrams of the type shown in Fig. |5l These are all diagrams of NRQCD with an ingoing 
QQ pair and a gluon and an outgoing QQ pair. They can involve singlet four-fermion 
operators and a gluon coupled to the quark or the antiquark line, but also four-fermion 
operators that couple to gluons. Such operators induce octet to singlet transitions on the 
QQ pair; they may be one of the operators listed in Eqs. ( [A.6| ) and ( [A.7 ), but also one 
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Figure 4: QCD Feynman diagrams describing the amplitude QQ g ^ 77 ^ QQ at leading order. 

of the four-fermion operators involving only covariant derivatives, which, despite being 
usually denoted as singlet operators, couple to the gluon field through the term —it"'gA°' 
in the covariant derivative. An example is the operator 

f 



OemC^o) = 3^^ 



-^D .ajx|0)(OlxM-^'D ■a)V', 



that annihilates (or creates) a.^Si color-octet pair and a gluon and creates (or annihilates) 
a ^Pq color-singlet pair through the term 



In the basis of Sec. 2.3.3, we obtain 



ImS8om(^5'o) = 0, 



Im/igjn(^5o) 



10 



a^QV, 



Im/iem(^S'o) = -a^Q'^vr, 
5 

Imtsemi^Pl) = 0, 
lmtscmCP2) = 0. 



(3.21) 
(3.22) 

(3.23) 

(3.24) 

(3.25) 
(3.26) 



The four-fermion operators to which the matching coefficients refer are listed in Appendix 
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Figure 5: Generic NRQCD four-fermion Feynman diagrams involving an ingoing QQ pair and a 
gluon and an outgoing QQ pair. The black box with a gluon attached to it and the empty box 
stand respectively for one of the four-fermion-one-gluon vertices and for one of the four-fermion 
vertices induced by the operators listed in Appendix ^ Eqs. (|A.3| )-( A.IO ). The black dot with a 
gluon attached to it stands for one of the quark-gluon vertices induced by the bilinear part of the 
NRQCD Lagrangian given in Eq. ( |A.l|) . 



Equations ( 3.2ip -( p^ ) are new. We see that the matching of the amplitude QQ g -^ 
77 allows to establish the values of the individual coefficients Im h em(^'S'o) and Im h cm(^'S'o)- 
Equations ( |3.24 ) and ( |3.26D disagree with those in Q. Equation ( |3.25| ) follows from the 
Landau- Yang theorem. 



3.4 QQg 



e ' e 



The matching of the QQg — > e~^e~ amplitude is performed by equating the sum of the 
imaginary parts of the QCD diagrams shown in Fig. ^ with the sum of all the NRQCD 
diagrams of the type shown in Fig. ^ 



In the basis of Sec. 2.3.3, we obtain 



ImS8em(^5'i) = -a'^Q^TT, 



Imi'U'Si 



1 

6 
23 



a^g^vr, 



54 



1 






ImK^i'^Su'Di 



1 

12 



a^Q^TT. 



(3.27) 
(3.28) 
(3.29) 
(3.30) 
(3.31) 
(3.32) 



The four-fermion operators to which the matching coefficients refer are listed in Appendix 



Equations ( 3.27D -( p!3^) are new. We see that the matching of the amplitude QQ g — > 
e^e~ allows to establish the values of the individual coefficients Im.h'g^{^Si), lm.h'^^(^Si), 
Im/igj^(^S'i,^ Di) and Imh'^^i^Si,'^ Di), whose sums were determined by the QQ —>■ e^e~ 
matching. 
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Figure 6: QCD Feynman diagrams describing the amplitude QQ g -^ e+e — + QQ at leading 
order. 

3.5 Cancellation of infrared singularities 

Once expanded in the gluon energy |g|, diagrams with a gluon attached to the external 
quark lines generate terms that are proportional to 1/|^. They come from the almost 
on-shell heavy-quark propagator that follows the gluon insertion. Such terms, which are 
singular for |g| — > 0, cancel in the matching. 

Let us consider two graphs like the first (second) and the fifth (sixth) in Fig. ^ or 
those in Fig. ^. Their sum is proportional to 

where e is the (transverse) polarization of the external gluon. Both quarks and gluon are 
on mass shell, therefore, we may write the denominator as 



{p + qY-NP 2{Ep\^-p-q)' 
By using the Dirac equation, we may also simplify the numerators 

v{p') ^ (-/ + M)= v(p")(/ + M) ^ - 2p' ■ ev{p') = -2p' ■ ev{p') 
[^ + M)^u{p) = ^{-^ + M)u{p) + 2p • £u{p) = 2p • £u{p). 

Then Eq. ( |3.33| ) becomes 

^ ^ 1 ^, ^ 1 



-p- ge 



En 



P-. 



. + p' ■ 9£ 



Ep'\q\ -p' -q 



1 + 0(1^-1) xTW 



(3.34) 



where M'^ is the amplitude calculated from the diagrams without gluon insertions but a 
color matrix t" inserted between the ingoing quark and antiquark. Equation ( ^.34 ) displays 
only terms that are singular for |^ ^ 0. 

In order to show the cancellation of the singular terms in the matching, it is sufficient 
to consider in NRQCD the last two diagrams of Fig. |5|, where the gluons couple to the 
quark and antiquark through 



-«.,!!|i!,a('i+ « 



M 



3 54 \ 

+ nT7J + --- V' + CC. 



2M2 8M4 



(3.35) 
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Equation ( [3.35 ) is the one-gluon part oi {D = d — it"' gA"") 



^^2-f = V'M:TirF + VJr + TinrFK + ---U + ^-«- (3-36) 






2M 8M3 16M5 



Vertices involving E or B fields are proportional to q. Hence, graphs with such vertices 
are not singular for \q\ -^ 0. In momentum space, the vertex induced by Eq. ( ^.35 ) on the 
quark line of Fig. g is 

«"^7l-^ + ^ + ...l=«"^1 (3.37) 



M \ 2M2 8M4 7 Ep 

similarly, the vertex induced on the antiquark line is 



-it"^—^. (3.38) 

Ep' 



The heavy-quark propagator in the second diagram of Fig. ^ reads 

^l + 0(|g1)); (3.39) 



Ep + \q\- Ep+g 1^ P-q 

Ep 



where Ep^g = \J M"^ + (p + g)^; the heavy antiquark propagator in the third diagram of 
Fig. ^ reads 

'l + 0(|g1)). (3.40) 



Ept + 1^ - -Ep'+g I ^ _ v' -q 



Ep, 



Finally, the sum of the two diagrams can be written as 

{-^'^' EM-p-i ^^'-' ^EM-f.^ (l + 0(k1)) xA<il„,cD, (3.41) 

where AI^j^qqj-, is the NRQCD amplitude calculated from the diagrams without gluon 
insertions but a color matrix t" inserted between the ingoing quark and antiquark. Equation 



(|3.4lD displays only terms that are singular for \q\ — > 0. 

Since the matching without external gluons guarantees that M. = A^nrqcd and also 



that M"" = -A^NRQCD' ^y matching Eq. (|3.34D with Eq. (|3.4l| ), all the displayed singular 



terms cancel. In general, the non singular terms do not cancel. Indeed, diagrams like the 
first (second) and the fifth (sixth) of Fig. ^ and those in Fig. y contribute to the matching 
coefficients of NRQCD. 

We have shown by an explicit calculation how infrared singularities cancel in the match- 
ing at leading order in a^ and to all orders in 1/M. The argument reflects the very general 
fact that QCD and NRQCD share the same infrared behaviour. 
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4. Conclusions 

We have performed the matching of the imaginary parts of four-fermion operators up 
to dimension 10. The matching enables us to write the electromagnetic S- and P-wave 
quarkonium decay widths into two photons or an e~^e~ pair up to order w^. In the power 



counting of Sec. 2.3.4, the result reads 



+^^^S^{HeS,)\QU'Si)\HeS,)) + ^^^{HeS,)\Q'U'S^)\HeS^)) 



27M6 ' ' -^/i-cmV -i;i V -i// g^6 

{HeS^)\V.„,eS^,'D,)\HeS^)) + ^-:^{hCS^)\%'ZCSi)\HCS^)), (4.1) 



■—^^^^{H{^Si)\'Pcm{^Si,^Di)\H{^Si)) + g^5 \l2y Oi)\^Q^^ 



TiHi'So) - 77) = ^^{HeSo)nmeSo)\HeSo)) 

-^-^^{Heso)\v.rneso)\Heso)) + ^^^^{Hesomu'someso)) 

+^^^{HeSo)\Q'U'So)\HeSo)), (4.2) 

FiHePo) - 77) = ^^^{HePo)n..ePo)\HePo)) 

'^"'^'"(i7(3Po)|^em(^Po)|^(^Po)) - ^^;P^{HePo)\TsU'Po)\HePo)), (4.3) 



5M4 



T{HeP2) - 77) = ^^^(i/(3P2)|Oem(^P2)|i^(^P2)) 



^^^^{HeP2)\VU'P2)\HeP2)). (4.4) 



Equation (4.1) has three extra terms with respect to the result of Ref. |^ (these are 
the terms proportional to {H{^Si)\Ssem{'^Si) \H{^Si)), {H{^Si)\Pen,CSi,^ Di)\H{^Si)) 
and {H{^Si)\ ^^^{^81) \H{^Si))), which are due to the different power counting. The 
other terms in Eq. (|4.1|) and Eq. ( [4.2[ ) agree with those in j^, where, however, only the 
sums of the matching coefficients in front of {H{^Si)\ Qem(^'S'i) iH^Si)) and {H{^Si)\ 
Q'U^Si) \H{^Si)), and {H{^So)\ Q'erai^So) \H{^So)) and {H{^So)\ Q'U^So) \H{^So)) 
were calculated. The coefficient of {H{^Pq)\ Vemi^Po) l-f^(^-fb)) iii (0) ^^ ^^^ same as 
in §, whereas the coefficients of (^(^Po)| Tsemi^Po) \H{^Po)) in ^) and of {H{^P2)\ 
Vemi^P2) \H{^P2)) in (4.4) are different. Equation (4.4) depends on two matrix elements 
less than the equivalent equation in j^. We could only partially trace back the origin of 
these differences. However, we stress that, at variance with the S-wave case, in the P-wave 
one, the differences with the previous literature cannot be reconciled by a different power 
counting. 

In the case of S-wave decays, Eqs. (|4.1| ) and ( |4.2| ) provide relativistic corrections with 
a few per cent accuracy in the charmonium case and with a few per mil accuracy in the 
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bottomonium one. At present, the largest uncertainties in the decay widths of the pseu- 
doscalars come from next-to-next-to-leading order corrections to the matching coefficient of 
the dimension 6 operator and next-to-leading order corrections to the matching coefficient 
of the dimension 8 singlet operator, which are both unknown. They are about 10% in the 
charmonium case and about 5% in the bottomonium one (in accordance to the values of 
as{M) and v listed in the introduction). In the charmonium case, the theoretical uncer- 
tainty is well below the experimental one. In the case of the decay widths of S-wave vectors 
in lepton pairs, next-to-next-to- leading order corrections to the matching coefficient of the 
dimension 6 operator are known and the dominant uncertainty comes from next-to-leading 
order corrections to the matching coefficient of the dimension 8 singlet operator. It is about 
10% in the charmonium case and about 2% in the bottomonium one. 



In the case of P-wave decays, Eqs. (4.3) and (|4.4|) combined with the next-to-leading 



order expressions of Im/em(^-Fb) and Im/em(^-P2) [0] provide a consistent determination 
of the electromagnetic widths with an estimated theoretical error that is about 15% in 
the charmonium case and about 5% in the bottomonium one. In |l^, |I[, it was observed 
that the inclusion of the next-to-leading order corrections to the matching coefficients alone 
(i.e. without relativistic corrections) brings the theoretical determination of the ratio of the 
charmonium ^Pq and ^P2 decay widths into two photons rather close to the experimental 
one. This supports the assumption that relativistic corrections follow, indeed, the non- 
relativistic power counting in v and are not anomalously large. The same observation, 
made for hadronic decay widths [|14| , p, urges a complete calculation of the relativistic 
corrections for P-wave hadronic decays of quarkonium, which is still missing. A first partial 
analysis is in ||^]. 

Equations (4.3) and ( [4.4| ) depend on three and two non-perturbative matrix elements 
respectively. Spin symmetry relates {H(-^Po)\Vcm{^Po)\H{^Po)) to {H{^P2)\ Vcmi^P2) 
\HeP2)): 



{HeP2)\VU^P2)\H{^P2)) = {HePo)\V,r^ePo)\HePo)) (l + 0{v)). (4.5) 

The matrix elements {H{^Po)\ Oemi^Po) \H{^Po)) and (^(^^2)! Oem(^i^2) [^(^^2)) satisfy 
an analogous relation, which, however, is not useful at relative order v'^. Therefore, Eqs. 
( |4.3| ) and ( [4. 4]) depend on four independent matrix elements. Since, at present, they are 
poorly known, in this work we refrain from phenomenological applications. We remark, 
however, that a further reduction in the number of non-perturbative parameters may be 
achieved in low energy EFTs like potential NRQCD (pNRQCD), where, in the strongly- 
coupled regime, the NRQCD matrix elements factorize in the quarkonium wave function 
and in few correlators of gluonic fields |15|. The completion of the P-wave hadronic de- 



cay width calculation at order u^ in NRQCD is the first step towards a complete analysis 
of both electromagnetic and hadronic relativistic corrections in the framework of the pN- 



RQCD factorization. Besides direct lattice evaluations of the matrix elements |16|, such an 
analysis has the potential to constrain the non-perturbative parameters enough to provide 
theoretical determinations matching the precision of the data. 
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A. Summary of NRQCD operators 

The two-fermion sector of the NRQCD Lagrangian relevant for the matching discussed in 
Sec. ^ is: 

r - l\(-n , D^ , S-gB , {D-gE) a-[-iDx,gE] , {&f , {D\a ■ gB} 

'-' " ^ 1^ ° ^ 2M ^ 2M ^ 8M2 8M2 + 8M3 + SM^ 

^ .{D\a.[-^Dx,gE]} + J-AD\{D.gE)} + -^]i^ 



32M4^ > L ,y jj i 32M4L >v u ji ' IQM5 j 

+ c.c, (A.l) 

where cr* are the Pauli matrices, iDo = ido — t°- gA^, iD = iV + t"- gA"-, [Dx,E] = D x 
E - E X D, E' = F'^ and B' = -eijkF^^ /2 (ei23 = 1). We have not displayed terms of 
order 1/M^ or smaller and matching coefficients of 0(as) or smaller. Equation ( |A.1| ) can 
be obtained by matching the QCD scattering amplitude of a quark on a static background 
gluon field along the lines of Ref. |^^. Note that in [ [17| ] the matching coefficient of the 
operator {D^,a ■ [—iDx,gE\]/M'^ is off by a factor 1/2. The right value may be found in 
Ref. |1|. 

The general structure of the four-fermion sector of the NRQCD Lagrangian projected 
on the QCD vacuum is: 



/-em _ V^ Cem ^^(n) / . r,N 

^4-f - 2^ A/f<in-4'-^4-fem- \^-^) 



Here, we list the operators relevant for the matching performed in Sec. y ordered by 
dimension {D = D — D). 



Operators of dimension 6 
Operators of dimension 8 



(A.3) 



Verai^S^) ^-^^ f-^"^) X|0) (0|xV + H.C., 

7'em('5i)=^V'Vx|0)(0|xV('-^'D') V + H.c, (A.4) 
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ou'Pi) =V'^ {-\^) xlo)(o|^^ {~\^) '^' 

am('Pi) =\i^^ {-\^ X a^ x|0)(0!x^ (-^^ X a) V, 

58cm('5o) =^V'^5^ • '?x|0)(0|xV + H.c, 
58cm('5i) =^V'^5^x|0)(0|x^'?^ + H.C.. 
Operators of dimension 9 

r8em('5o) =^V'^xlO)(0|x^('^ -gE + gE- "d )V' + H.c, 
'^'selC'^i) =^V'^^X|0) (0!x^a^('^ -gE + gE. V)i; + H.c, 
rg^emC'^i) =^V'^^x|0)(0|x^(T X (-"d xgE-gEx "d)^ + H.c, 
'?'8em('5i) =\^^dx\^){^\x^^ x(d xgE-gEy. "d)^ + H.c, 
'?'8e'm(''5i)=^V^^^x|0)(0|xV^-(^M^-)+5^(^'^^'))V + H.c, 

T8cm('Po) =^V'^ p • d) X|0)(0|X^^ • ^^V' + H.C, 

T8em('Pi) =^V'^ p X a) x|0)(0|x^<? X 5^^ + H.c, 

T8em('P2) =^V'^ (^D^V^)) xlO) (01x^(^5^5;^^ + H.C. 

Operators of dimension 10 

Q"cm('5o) =^^t ("i"^) x|0)(0|xV + H.c, 



Qcm(''5i) =\^^ ("i"^) <5^xlO)(0|x^aV' + H.c, 






(A.5) 



(A.6) 



(A.7) 



(A.^ 
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ru'Pi) =li^^ (-^'^) {-'2^') ^lo)<oi^^ (-^^^) v + H.c, 

VomCPo) =^V'^ (-'^^D ■ A (-'^^) X\0){0\X^ (-'^^ ■cj]i; + H.C., 



^em(=^P2) =^V'^ (-0 '^(V-') (-^'^) X|0)(0IX^- (-^"^(V)) V + H.C, 






B. Summary of matching coefRcients 

In the following, we list all imaginary parts of the matching coefficients of the electromag- 
netic four-fermion operators up to dimension 10, calculated at 0(1) in the strong coupling 
constant in Sec. |^. The matching coefficients refer to a basis of operators from where 
T8{^So), Tg (^S'l), Tg {^Si) and Tg (^Si) have been removed by suitable field redefini- 
tions. 



Operator of dim. 6 


Matching coefficient 


Im (Value) 


O^mCSo) 


Imf^^CSo) 


a^QV 


0,mCSi) 


Im/en.('5l) 


■'- 2^2 

3a Q^ 
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Operator of dim. 8 

^cm('5i) 



^cm( -S*!, Di 



Matching coefficient 

Im5rcm(^5'o) 
Im5rcm(^<S'i) 



Im 5cm (^ 5*1,^1)1 



Im (Value) 



a^Q\ 



a^oV 



ra^QV 



Ccm( -Pi) 
Ccm( -Po) 
Ccm( -Pi) 

Oem('i^2) 



Im/e,n('Pl) 
Im/e,n('Po) 
Im/e,n('Pl) 

Im/e„,(3p2) 








'58cm( 'S'o) 
'58cm(' 'S'l) 



ImS8cm(^5'o) 
ImS8cm(^5'i) 



-a^Q\ 



Operator of dim. 9 

-'SemV '-'IJ 

^cm( Pb) 

^cm( P'l) 
'?8cm( P2) 



Matching coefficient 

Im41(35i) 

Imt8cm(^Po) 

Imt8cm(^Pl) 
Imt8cm(^P2) 



Im (Value) 



rO^OV 



:a^Q\ 








Operator 


of dim. 10 


Matching coefficient 


Im (Value) 


QcmC^o) 




Irah'U^S^) 




10 2r,4 


QemC^o) 




luih'U^S^) 




¥'<}'' 


Qem('5l) 




lrnh',^{^Si) 




i"^Q^- 


QemC^l) 




Irnh'U'Si) 




ia.^g^. 


Q^mC^l,^ 


^1) 


Irnh'U^S^? 


Di) 


r''}'' 


QemCSl,^ 


^1) 


Irnh'U^S^? 


Di) 


h-"^'^ 
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Von^i^Pl) Im^emCPl) 

Vemi^Po) Im5em('Po) -7a^Q\ 

V,rn{^Pl) Im^emfPl) 

5 
20 

Vem{^P2?F2) luiQ.^i^ P2? F2) "^"'^^^ 

15 

Qcm('^I)2) Im/lcm('^I)2) 

Qcm(^I)3) Im/lem(^I)3) 
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